Singularities in ground state fidelity and quantum phase transitions for the Kitaev model 
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The ground state fidelity per lattice site is shown to be able to detect quantum phase transitions for the 
Kitaev model on the honeycomb lattice, a prototypical example of quantum lattice systems with topological 
order, h is found that, in the thermodynamic limit, the ground state fidelity per lattice site is non-analytic at 
the phase boundaries: the second-order derivative of its logarithmic function with respect to a control parameter 
describing the interaction between neighboring spins is logarithmically divergent. A finite size scaling analysis 
is performed, which allows us to extract the correlation length critical exponent from the scaling behaviors of 
the ground state fidelity per lattice site. 
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With the advent of its discovery in the fractional quantum 
Hall effect [1], topological order emerges as a new paradigm 
in the study of quantum phase transitions (QPTs) |2]. Sub- 
sequent investigations show that topological order occurs in 
various strongly correlated lattice systems undergoing QPTs. 
A characteristic feature of quantum systems with topological 
order is their insensitivity to any local perturbations [3 J. Such 
an essential difference between topological and symmetry- 
breaking orders invalidates the usual tools used to describe 
a symmetry-breaking order, such as long range correlations, 
broken symmetries, and local order parameters 13] • 

Recently, much attention has been paid to an exactly solv- 
able spin 1/2 model on a honeycomb lattice introduced by 
Kitaev Jstl for fault-tolerant topological quantum computa- 
tion 10]. The model describes a set of spins located at the 
vertices of a two-dimensional honeycomb lattice, subject to a 
spatially anisotropic interaction between neighboring spins. It 
has been shown that it carries excitations with both Abelian 
and non-Abelian braiding statistics, which do not obey or- 
dinary bosonic and fermionic statistics, but are anyons with 
more intricate statistical behavior [7]. An experimentally fea- 
sible realization of the model in a system of cold atoms on an 
optical lattice has been addressed |8(] (see also with 
the expectation to perform quantum computation by utilizing 
braiding of collective excitations implanted in topologically 
ordered coherent quantum many-body states. 

On the other hand, a viable scheme to determine the 
ground state phase diagram of a quantum lattice system with- 
out prior knowledge of order parameters was proposed in 
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14tl . This was achieved by studying the 
singularities in the ground state fidelity per lattice site ] 15]. In 
fact, the ground state fidelity may be interpreted as the par- 
tition function of a classical statistical vertex model with the 
same lattice geometry by using the tensor network representa- 
tions of quantum many -body wave functions ] 12]. Therefore, 
the fidelity per lattice site is nothing but the partition function 
per site in the classical statistical vertex lattice model 11^. 
This justifies why QPTs may be detected as singularities in 
the fidelity per lattice site as a function of the control param- 
eters (see also Refs. 11171 llSll for the connection between the 



fidelity and QPTs). Therefore, an intriguing question is to see 
if the fidelity approach captures the physics underlying QPTs 
in quantum lattice systems with topological order 

The purpose of this paper is to show that the ground state 
fidelity per lattice site is able to detect QPTs for the Kitaev 
model on the honeycomb lattice, a prototypical example of 
quantum lattice systems with topological order. First, we 
derive the ground state fidelity per lattice site between dif- 
ferent ground states from the exact solution of the Kitaev 
model on the honeycomb lattice. This is achieved by ex- 
ploiting the fact that the original spin model on the honey- 
comb lattice is rephrased as a p-wave BCS model with a 
site-dependent chemical potential for spinless fermions on a 
square lattice d (see also Refs. iH IH). The ground 
state of the latter is a BCS type state, as a consequence of 
the Jordan- Wigner, Fourier and Bogoliubov transformations. 
Second, the phase boundaries separating the gapless phase 
from different gapful phases are reproduced by investigating 
the singularities in the fidelity per lattice site as a function 
of the control parameters. It is found that, in the thermody- 
namic limit, the ground state fidelity per site is non-analytic 
at the phase boundaries. That is, the second-order derivative 
of its logarithmic function with respect to a given control pa- 
rameter is logarithmically divergent as the phase boundaries 
are crossed. Third, we perform a finite size scaling anal- 
ysis for the Kitaev model, aiming at extracting the correla- 
tion length critical exponent from the scaling behaviors of the 
ground state fidelity per site. Our exact results offer a bench- 
mark to investigate QPTs for two-dimensional quantum lattice 
systems with topological order numerically in the context of 
tensor network representations 123^ ,24. ,25. ,26J . 

The Kitaev model on a honeycomb lattice. Consider a spin 
1/2 model on a honeycomb lattice with the Hamiltonian ]|5|] 
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x- bonds 



V- bonds 
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-bonds 



where Ja are interaction (control) parameters and cr". are the 
Pauli matrices at the site j, with a - x,y and z- The Hamil- 
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Ionian ([T]l may be fermionized by performing the Jordan- 
Wigner transformation |T3 2^ 22 1 from the PauH spin 
matrices cr" to the spinless fermion operators c\ and cj. This 
one-dimensional fermionization is realized by deforming the 
hexagonal lattice into a brick-wall lattice which is topolog- 
ically equivalent to the original lattice. If we introduce the 
Majorana fermions: A = (c - c' )// and A - c + c' , then the 
Hamiltonian ^ becomes 

H = -/y, ^ A,,Ah + iJy ^ AiA„. 

x-bonds y-bonds 

-iJ, ^ arAbA„, (2) 

z-bonds 

where the subscripts w and b denote two sublattices in the 
brick-wall lattice, and = Bi,B„ along the z-bond is con- 
served lEoil . with r being the coordinate of the midpoint of the 
bond connecting the /j-type and w-type sites. This in turn is 
equivalent to a model of spinless fermions on a square lattice 
with a site-dependent chemical potential: 

H = J,Yj{4+d,){dl.-dl,) 

r 

+j^.Y^(di+d,)(di^-di,) 

r 

+J,Y^ar(2dU,.-l). (3) 

r 

Here the unit vector and e, connects two z bonds and 
crosses a x- and y-bond, respectively. For large enough sys- 
tems, the ground state configurations are bulk vortex-free con- 
figurations isl I21I], which implies a,- - 1 for all r. There- 
fore, the ground state may be obtained by performing a fourier 
transformation. Up to an unimportant additive constant, the 
Hamiltonian (|3]l in the vortex-free configuration now reads, 

Hg = (f*^'^* + Y i4d'lk + H-C.)J , (4) 

with €k = 2J, - 27 V cos A:^ - 27,, cos ky, and = 27,. sin q^ + 
2JySmky. The Hamiltonian (|4]l is a p-wave type BCS pairing 
model and can be diagonalized by means of the Bogoliubov 
transformation. It yields that the BCS type ground state is 
l^> = Ukiuk + Vkdldi^m, where = 1/2(1 + eulEk) and 
\vk^ - 1 /2(1 - eitlEk), with the quasiparticle excitation energy 

Ek = ^e2 + A2|ll9fi. 

The ground state fidelity per lattice site. Consider two 
ground states \g) and \g') corresponding to different values of 
the control parameters 7 = (7v, 7y, 7^) and 7' = (7^, 7,', 7!), re- 
spectively. The fidelity F{ J; J') = {g'\g) asymptotically scales 

as F{ J; J') ~ d(J; J') , with the total number of sites in the 
lattice. Here d{J; J') is the ground state fidelity per lattice site, 
introduced in Refs. 1 11, 12]. Although F{J; J') becomes triv- 
ially zero for continuous QPTs, the fidelity per lattice site is 
well defined in the thermodynamic limit: 



It satisfies the properties inherited from the fidelity F(J; J'): 
(i) normalization d(J; f) - I; (ii) symmetry d{J; J') - d(J'; 7; 
and (iii) range < d{f; /')<!. 

For the Kitaev model on the honeycomb lattice, the loga- 
rithmic function of the fidelity per site. In dh{J', 7'), is half of 
the logarithmic function of the fidelity per site, In d^qi J; J'), 
for the model of spinless fermions on a square lattice. This re- 
sults from the fact that the number of sites in the honeycomb 
lattice doubles that of sites in the square lattice. The BCS type 
ground state \g) yields the ground state fidelity per lattice site 
for the spinless fermion model on the square lattice: 

In dsqi J ; f) ^ -T—:^ dk, dk,.ln(ulu[ + vWA, (6) 
(2ny Jo Jo 

where Uk and Vk depend on 7, whereas and depend on 
J'. Here we emphasize that although the information about 
the topological nature of the Kitaev model is lost in the spin- 
less fermion representation, the unitary equivalence between 
the two representations preserves the fidelity. Since the extra 
prefactor does not affect the singularities in In dh(J; J') and 
\ndsqiJ; J'), hereafter we focus on \ndsqiJ; J') to carry out 
the scaling analysis below, and omit the subscripts for brevity. 

For a finite-size system, the Hamiltonian (|4]i, resulted from 
the Jordan-Wigner, Fourier and Bogoliubov transformations, 
depends on boundary conditions imposed on the original spin 
model ([TJ. In contrast to open boundary conditions, there is an 
extra boundary term if one adopts the periodic boundary con- 
ditions. However, such a boundary term does not contribute 
to the fidelity per site, although it carries the topological de- 
pendence of the ground state degeneracy 1 19]. From now on, 
we are only concerned with the fermion model on a square 
lattice with the periodic boundary conditions (i.e., a torus) to 
analyze the ground state fidelity per lattice site for finite-size 
systems Ii27jl . from which it is sufficient to extract the bulk be- 
haviors of the model. As such, for a system on a torus with 
an even linear size L, the logarithmic function of the ground 
state fidelity per lattice site, In d( J; J'), takes the form: 



\nd(J- J^)^j2Yj ^^"'k^-'k + "l^'k)- 

k„k. 



(7) 



diJ~-f) = lim FHJIJ')- 

N—><x> 



(5) 



Here kx and ky take values from the set: 7Tm/L(m - -(L - 
l)/2, . . . , (L- 1)/2), and the double summation is over all pos- 
itive values of both k^ and ky. 

Ground state phase diagram and singularities in the ground 
state fidelity per lattice site. Now we turn to the ground 
state |)hase diagram. This follows from the singularities in 
\nd{J,J'). One may show that lnii(7;7') in Eq. ^ and 
the first-order derivative with respect to a control parame- 
ter is continuous, but the second-order derivative logarith- 
mically diverges when the phase boundaries determined by 
17,1 = 17,1 + 17,1, 17,] = 17,1 + \JA and |7,| = |7J + |7,] are 
crossed. This is consistent with the original analysis by Ki- 
taev 15] (see also Refs. | liSIHl]). In Fig.fTIa), we plot the 
logarithm of the fidelity per site, In d{J; J'), as a function of 



7r and J[ for Jy - J^- 1/2 and J[, - Jl - \ 12. It is seen that 
a pinch point occurs at {JxcJxc) - (1, !)■ That is, there are 
singularities along the lines - 1 and J'^- 1 . Therefore, the 
drastic change of the ground state many-body wave functions 
at 7vc is reflected as the singularities in lnii(7; J'). Similarly, 
the numerical results are plotted in Fig.lTfb) for the logarithm 
of the fidelity per lattice site. In d{J, J'), as a function of /, 
and J[ for fixed Jx - J'^ - Jx - J'y - 1/2, with a pinch point 
ata„7^,) = (l,l). 




FIG. 1 : (color online) (a) The logarithm of the fidelity per lattice site, 
In d(J, J'), is shown as a function of 7,- and y(. for fixed Jy = i(, = 
J- = J'. = 1/2. It exhibits a pinch point at {JxcJxc) = (\, !)■ (b) 
The logarithm of the fidelity per lattice site, \ad( J , J'), is shown as a 
function of J- and J'- for fixed 7, = J[ = Jy = Jy = 1/2. It exhibits 
a pinch point at {J-c, J^c) = (1> 1). Here a pinch point is defined as an 
intersection of two singular lines. 

More precisely, for any fixed J', In d(J; J') is logarithmi- 
cally divergent when J' is varied such that a critical point is 
crossed. Suppose Jy and J^ are fixed, and only Jx is a control 
parameter that varies. Then we have 



dHndjJ, J') 
dJ? 



= ki In 17,: - Jxc\ + constant, 



(8) 



where ki is a non-universal prefactor that depends on Jy, J^ 
and J', and Jxc is the critical value of 7, for fixed 7, and 7-. 
The numerical results are plotted in Fig. |2a) for 7,. = 7,- = 1 /2 
and Jxc - 1. The least square fit yields ki 0.02360. Simi- 
larly, we have presented numerics in Fig. |2tb) for the second- 
order derivative of In d(J; J') with respect to 7,, with Ji - 0.8 



and Jx = Jx 



J'y 



1/2. It turns out that it diverges 



logarithmically in the same way as (O with Jx replaced by 7,-, 
and ki « 0.04726. 

Finite size scaling analysis. For a system of finite size 
N = L? (with L the linear size), there is no divergence in 
the second-order derivative of In d{J, J') with respect to Jx, 
since QPTs only occur in the thermodynamic limit. Instead, as 
seen in Fig. I2a), some pronounced dips occur at the so-called 
quasi-critical points 7„„, with the dips values logarithmically 
diverging with increasing linear size L, 



dHnd(J, 7') I 
dJ? I 



A;2 In L -H constant. 



(9) 



where k2 is a non-universal prefactor k2, which takes the value 
k2 ~ -0.02312 for 7; = 0.8 and 7,, = J'y = 7, = 7; = 1/2 




FIG. 2: (color online) (a) The second-order derivative of the log- 
arithm of the fidelity per lattice site, lnd(J,J'), with respect to 7, 
diverges at the critical point in the thermodynamic limit. However, it 
remains analytic for finite-size systems, although more pronounced 
dips occur with increasing linear system size. Here J[ = 0.8 and 
7,, = 7,, = J^ = J'. = 1/2. (b) The second-order derivative of the 

logarithm of the fidelity per lattice site, In d{J, J'), with respect to 7; 
diverges at the critical point in the thermodynamic limit. However, it 
remains analytic for finite-size systems, although more pronounced 
dips occur with increasing linear system size. Here 7' = 0.8 and 
7,- = J[ = Jy = J'y = \ 12. (c) The dips values scale as In L with 
increasing linear size L for 7^ = 0.8 and Jy = 7', = J- = J'. = \ jl. (d) 
The dips values scales as InL with the linear size L for 7' = 0.8 and 
7, = 7J. = 7v = 7,'. = 1 12. (e) The positions of the dips approach the 

critical point Jxc = 1 with increasing linear size L. Here d(J, J') is 
shown as a function of Jx for 7J. = 0.8 and 7,, = J'y = J, = J'. = 1/2. 
(f) The positions of the dips approach the critical point J-c = 1 with 
increasing linear size L. Here d(J, J') is shown as a function of 7; for 
7-' = 0.8 and 7,- = 7; = Jy = J'y = 1/2. 



(see Fig.|2lc)). In addition, Jxm approaches the critical value 
as Jx„ ~ 1 - 3.96384L-i '"^245^ follows from Fig.|2e). 
The scaling ansatz in the system exhibiting logarithmic di- 
vergences requires that the absolute value of the ratio ki/k2 
is the correlation length critical exponent v. In this case, 
\ki/k2\ ~ 1.02076, very close to the exact value v - I. This 
is consistent with the fact that the gap A for the Bogoliubov 
quasiparticle scales as A ~ 7^; - Jxc near the critical point 
Jxc- Similarly, a finite size scaling analysis is performed for 
In d{J, J') with 7^ = 0.8 and Jx = J'y = Jy ^ J'y ^ 1/2. 
In Fig. lite), flie least square fit yields k2 « -0.04640. The 
numerics for d^lnd(J,J')/dJx^\j,=j,,„ and 73,, are plotted in 
Figs.|2|d) and (f), respectively. 

In order to address the scaUng ansatz for a system exhibit- 
ing logarithmic divergence fl^, we take into account the dis- 
tance of the minimum of In d(J, J') from the critical point 
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to investigate 1 - exp[5y lnd{J, J') - lnd{J, /')|y,=y„J as 
a function of L(J^ - J^m) for different linear sizes Us. The 
numerical results for the linear size ranging from L = 401 up 
to L = 1401 ai-e plotted in Fig.^a). All the data for differ- 
ent L's collapse onto a single curve, indicating that the model 
is scale invariant, i.e., %IL - "B,' jL' , and that the correlation 
length critical exponent v = 1. The same conclusion can be 
drawn from Fig.|3lb), where the data collapsing is confirmed 
for 1 - exp[52^ In d{J, J') - d\ In d( J, /')|y,=y,J ■ 
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FIG. 3: (color online) (a) A finite size scaling analysis is per- 
formed for a quantity defined as D{J, J') = \ — exp[35 \nd{J, J') - 
d] ln<i(7, 7')Ia=j.™]' Jy,J- and J' fixed. The scaling ansatz 
for logarithmic divergences implies that D{J, J') is a function of 
LiJy - for fixed Jy, J- and J' . (b) A finite size scaling analysis is 
performed for a quantity defined as D{J, J') = l-exp[3y lnrf(7, J')- 
d'j \nd{J, y')ly-=A„,]! with J^, Jy and J' fixed. The scaling ansatz im- 
plies that D{J, J') is a function of L{J. - J,^) for fixed /j, Jy and J' . 
All the data from L = 401 up to L = 1401 collapse onto a single 
curve. This shows that the system at a critical point is scale invariant 
and that the correlation length critical exponent v is 1 . 

Summary. We have demonstrated that the ground state fi- 
deUty per lattice site is able to detect QPTs in the Kitaev model 
on the honeycomb lattice. It is found that, in the thermody- 
namic limit, the ground state fidelity per lattice site is non- 
analytic at a critical point. More precisely, the second-order 
derivative of its logarithmic function with respect to a given 
control parameter is logarithmically divergent as the phase 
boundaries are crossed. A finite size scaling analysis has also 
been performed to extract the correlation length critical expo- 
nent from the scaling behaviors of the fidelity per site. Our ex- 
act results offer a benchmark to numerically investigate QPTs 
for two-dimensional quantum lattice systems with topological 
order in the context of tensor network representations, which 
is currently under investigation. 
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